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1. Introduction 



Radiative transitions play an important role for our understanding of QCD, in particular of 
heavy quarkonia. They provide information about the wave functions describing the physical sys- 
tem and probe both the perturbative and non-perturbative regime. Especially El transitions give 
significant contributions to the total decay rate and yield clean signals, which are observed in the 
experimental facilities. In the last few years CLEO, BES and the B factories have improved their 
observations of radiative transitions, a review about recent developments can be found in M\. 

On the theory side, electric dipole transitions were treated in several potential models, a sum- 
mary can be found in We will refer to [Q] for comparison with our results. A model-independent 
treatment to check and improve the calculations has been missing so far. However, in the last decade 
there has been significant progress using effective field theories (EFTs) to describe heavy quarko- 
nium (see and references therein). Since heavy quarkonium is assumed to be a non-relativistic 
system we may take advantage of the hierarchy of scales m^>mv^$> mv , where v <C 1 is the heavy 
quark velocity, m is the heavy quark mass ("hard scale"), p ~ mv is the relative momentum of the 
bound state ("soft scale") and E ~ mv 2 is the binding energy ("ultrasoft scale"). The ultimate EFT 
living at the ultrasoft scale is potential non-relativistic QCD (pNRQCD). In 2005, for the first time 
radiative decays, concretely Ml transitions, were calculated in this theory Using the frame- 
work of that paper as a guideline we close the remaining gap and compute the decay rates of El 
processes between S and P states (like n 3 Pj -> n'^Si y). The following is based on 



2. Framework 



By integrating out the hard scale m S> Aqcd from the fundamental theory (QCD) in perturba- 
tion theory (a s (m) < 1) one obtains non-relativistic QCD (NRQCD) [[7], 0]. For the calculation of 
El transitions at relative order v 2 only the two-fermion Lagrangian S£i-f matters and the relevant 
part reads 

< D ° + 2Tn + 5?) V + Ce ^ (£ CT - B "" +i S? ff -[ DX ' E ""l) " + " - <2 " 

with /D = id - gT a A a Q - ee Q A e ^, /D = z'V + gT a A a + ee Q A em and y denoting a Pauli spinor for 
the heavy quark. The matching coefficients are found to be eg" = 1 +C F a s {^ H )/1n + 0{c%) and 

c em = 2c em _ j wkh ^ ^ m 

For processes at the ultrasoft scale, NRQCD is not yet the appropriate theory, since there 
are still several scales entangled (p,E,Aqcd) and thus no homogeneous power counting can be 
established. Integrating out the soft scale mv we obtain a theory for ultrasoft modes, i.e. pNRQCD 



[q- llPJ- The crucial step to disentangle the energy and momentum scale is the multipole expansion 
in the relative distance r. To be definite we will use the power counting in the weak-coupling 
regime (p> E > Aqcd), which reads 

r ~ 1/jbv, V r = d/dr~ mv, V = d/dR ~ mv 2 , E,B ~ (mv 2 ) 2 , E em ,B em ~ k 2 . (2.2) 

ky is the energy of the emitted photon, which scales like mv for transitions between states with 
different principal quantum numbers. 
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The pNRQCD-Lagrangian contributing at NLO in the decay rate, i.e. at order kyV° /ni 2 , reads 



+ gV A (O t r-E5 + 5 t r-EO)} 

+ -^ypNRQCD + ^ ight , (2 - 3 ) 

where the covariant derivatives are given by iD§0 = idoO—g[T"AQ, O] and /DO = iVO+g[T a A a , O] 
and the trace goes over the color and spin indices. The singlet potential Vs has been calculated per- 
turbatively and non-perturbatively to order 1 /m 2 (| 11, 12, 13Q , for more original references see [0), 



we display the structure of the relevant potentials for computations at NLO in the decay rate, 

Ks(r)= y<»> (r) + ^ + !fM + !%>, (2 . 4) 
m m z m z 



VS to = V S WL • S + 4 2) (r)S 2 + y£> (r) [3(f • 0l ) (f • <7 2 ) - ffi ■ <T 2 ] . (2.6) 



Vi 2) (r) =V r (2) W+ rfW,? 2 }!^, (2.5) 



The relevant part of Jz^nrqcd for El transitions is 



-^ypNRQCD = ee Q J d 3 rTr^V rE Sh-E em S + VS E Oh-E em O + ^V^ rV)2rE S' f r- [{rV) 2 E em ]S 



+ jJ_yV-(rxB) 5 tr V rxB em }S 
Am 



+ i -LyM<-x(<-v)fi) 5 t{v r . ; r x [(rV)B em ]}5 



12m 

+ J_y('V)t7-Br 5 t CT i . [( rV )B em ]5 
4m 

(£XV,) [5t ' a] ' (r '" X Vr)5 } ' (2 ' 7) 

In fact more terms are allowed according to the symmetries of pNRQCD. However, we can show 
that their matching coefficients vanish. The matching is done by equating Green's functions in 
NRQCD and pNRQCD at the energy scale mv order by order in the inverse mass. 

The crucial argument for several operators is that diagrams in NRQCD which can be cast into a 
reducible structure also give reducible diagrams in pNRQCD . Therefore they have to be subtracted 
to obtain irreducible operators in pNRQCD and do not play a role in the matching procedure. An 
example is the diagram in Fig. 1, where the gluonic contribution can be factorized out yielding just 



a potential. Using this argument we can fix all of the Wilson coefficients in (2.7) at leading order 



in a em , so that the exact results reproduce the ones from tree level calculations, namely 

yr-E _ yr-E = y(rVfr-E = yV-(rxB) = yV r -(rx(rV)B) = ^ y(rV)a-B = £ em ^ yO-(ExV r ) = Q em _ ^.8) 

With the help of the formalism developed in [Q| we can describe the states in a quantum 
mechanical way using wave functions and compute the decay rate at NLO, i.e. at relative order v 2 , 
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Figure 1 : Example for a reducible diagram, if the electromagnetic operator commutes with the gluonic ones. 
It does not contribute to the matching coefficient of a single operator. 



from the Lagrangian (2.7). We obtain 



where 



/(/+!] 



(2.9) 



IN 



poo 

/ drr»R n , Q {r)R nl {r). (2.10) 
Jo 



R contains all of the wave-function corrections due to the higher-order potentials mentioned in 
(2.4)-(2.6), the relativistic correction of the kinetic energy, — p 4 /4m 3 , and higher-order Fock state 



contributions due to intermediate color-octet states. In contrast to Ml transitions the latter ones do 
not vanish for El decays (see [||] for explicit expressions). 

The expression ( |2.9[ ) is also valid in the strongly coupled regime (without color-octet contri- 
butions in R), where p ~ Aqcd, since we made use of non-perturbative matching arguments and 
additional operators do not appear in this regime. 

Compared to the results with the potential model calculation in [^] we find an equivalence 
between ( |2~^ ) and the corresponding formula there at the given order. However, our definite power 
counting allowed us to include all relativistic corrections systematically, in particular the color-octet 
contributions in the weak-coupling regime and the one coming from the potential Vr l \ Both were 
missing in former approaches. Furthermore we can show that the anomalous magnetic moment 
c e p n — 1 ~ 0(a s {m)) is actually suppressed and does not lead to large non-perturbative contribu- 
tions. 

Without much effort one can extend the discussion to other processes like n 1 P l -> n' x S y 
and n 3 S\ — > n' 3 Pjy (see also for transitions between states with the same principal quantum 
number, where corrections ~ k r are suppressed. 



3. Phenomenological analysis 

Based on these results a phenomenological analysis for bottomonium and charmonium de- 
cays at relative order v 2 can be performed. For a complete analysis we need a parametrization of 
chromoelectric field correlators in the weak coupling regime. Since El transitions always involve 
excited states, we alse require the quarkonium potentials in the strong coupling regime. These have 
to be matched with the known short distance behaviour. 
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As a first approach we proceed as following: As a parametrization of the static potential at 
short distances we use the perturbative expression at 3 loop with leading ultrasoft resummation 
with the parameters given in jl4| , p~5| ] (see also references therein) 1 derived from a matching of 
the static energy to unquenched lattice data [|l7|]. For large distances we use the string potential 



Vstring = —7t/l2r + ar + C [ |18| ] matched to lattice simulations [ ]17[ ] at r = 1.5ro (where r$ is the 
Sommer scale). These two potentials merge together smoothly at r ~ 0.8ro. To obtain the leading 
order wave functions we solve the Schrodinger equation with the static potential numerically using 



the Mathematica program schrodinger.m [19]. We fix the charm and bottom mass in our scheme a 
posteriori by matching the obtained masses for the J ' j\\f and T(15) with the physical ones. 

Concerning relativistic corrections the main contribution arises from wave function corrections 
due to subleading potentials. For short distances, here for r < (2GeV) _1 , we apply perturbative 
results at LL, wheras for long distances, r > (2GeV) _1 , we use parametrizations from quenched 
lattice results [20|, 21, [22| ] as a non-perturbative input. 2 Future approaches should aim for smooth 
transitions between these two regimes. We neglect color octet effects, which cannot be determined 
from current lattice simulations. 

The results of this computation are given in table [j] for bottomonium and in table ^| for char- 
monium decays. We do not include decays with ky > (p), where our power counting is assumed 
to break down. We see that the relativistic corrections lower the decay rates considerably, by 10- 
30% for bottomonium and by 20-60% for charmonium, which is especially striking for the decays 
h c (lP) ->■ ri c (lS)j and y^(2>S) — > % c o(lP)Y- The wave function correction due to the potential Vr 
yields particularly large contributions. The expansion works much better for bottomonium, since 
the average relative velocity is smaller {v\ ~ 0.1, v 2 c ~ 0.3). We estimate the uncertainty of our 
NLO result to be of order 10% for bottomonium and of order 30% for charmonium. This is on 
the one hand the generic size of a correction at 0(v 2 ), which can arise from color-octet effects or 
a systematic error in the treatment of the subleading potentials. On the other hand this is also a 
conservative measure for the total higher order effects, which are supposed to be suppressed by v 
compared to the total NLO corrections. 

Comparing our values with the potential model calculations in [Q] (for scalar and vector con- 
fining potential) we tend to get slightly larger values for bottomonium and slightly smaller values 
for charmonium. Within our uncertainties we stay consistent with observations, provided both the 
branching fraction and the total decay rate have been measured. 

We emphasize that except for the masses of the J / \j/ and Y(1S) and the photon energies no ex- 
perimental input has been used. Additional input and a thorough investigation of the fine splitting 
effects in the spectrum would improve the predictive power of our results. This analysis should 
be repeated, when unquenched results for the required chromoelectric field correlators and espe- 
cially for the subleading potentials become available, with emphasis on a proper connection to 
perturbative results, maybe at higher order (NLL), and a more elaborate uncertainty estimate. 



'To obtain a convergent behaviour one has to perform a renormalon subtraction at a scale p [jlfj]. In Jh] ] a s (l/r) 
was expanded in terms of oc s (p) to get a more stable behaviour. For short distances, r < 0.14ro, large logarithms yield 
an unreasonable shape of the potential, therefore we expand a,(p) in terms of a s (l/r) in this regime. 

2 As far as available we use fits with parametrizations based on calculations from the string model [H]. ( so f ar 
(2) — 
undetermined) and V^ 2 (small) are set to in the non-perturbative regime. 
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process 


r pNRQCD/ keV 






ltx D p G /keV 


AbO\ Lr J 




31.8 


29.7 ±3.1 


25.7-27.0 


- 


Yh-i (\P) 
Abl\ Lr J 




40.3 


35.8 ± 4.0 


29.8-31.2 


- 


Yu~> ( 1 Pi 


— \ T(]S)v 


45.9 


40.6 ± 4.6 


33.0-34.2 


- 


hu(]P) - 


-s- ni (\S)v 


60.8 


44.3 ± 6.1 


- 


- 


YOS,\ — $ 

1 \/-KJ J 7 


y, „ ( 1 P") v 

X,M)\>- r ) / 


1.52 


1.13 ±0.15 


0.72-0.73 


1.22 ±0.16 


T(2S) -) 




2.26 


1.94 ± 0.23 


1.62-1.65 


2.21 ± 0.22 


Y(2S) -» 


Xbi(iP)y 


2.34 


2.19 ±0.23 


1.84-1.93 


2.29 ± 0.22 


Xbo(2P) 


^T(25) 7 


12.6 


13.0 ± 1.3 


10.6-11.4 




Xbi(2P) 


^T(25) 7 


17.1 


16.3 ± 1.7 


11.9-12.5 




Xbi{2P) 


^T(25) 7 


20.4 


18.1 ±2.0 


12.9-13.1 




T(35) -» 


Xbo{2P)y 


1.44 


1.05 ±0.14 


1.07-1.09 


1.20 ±0.16 


T(35) -» 


Xbi(2P)y 


2.38 


2.05 ± 0.24 


2.15-2.24 


2.56 ± 0.34 


T(35) -» 


Xbi{2P)Y 


2.53 


2.35 ± 0.25 


2.29-2.44 


2.66 ± 0.41 



Table 1: El decay rates for bottomonium. Our pNRQCD results compared to a potential model calculation 
[[3J and and the current PDG values [p4|]. LO denotes the result obtained without relativistic corrections, 
NLO indicates the result up to ff(v 2 ) neglecting color-octet effects in the weak-coupling regime and non- 
perturbative contributions to v/ 2 '. The error estimates give the generic size of one ff{\' 2 ) correction as well 
as an estimate for the sum of all corrections at <^(v 3 ). For hy{\P) — > ri b (lS)Y we have taken '^(is) = 9402 
GeV p5| ] to determine the photon energy. 



process 


1 pNRQCD v 


r NLO /keV 
1 pNRQCD v 


r! d /keV 


r™ G /kev 


Xco(lP) 




199 


158 ± 60 


162-183 


122 ± 11 


XdM 




421 


302 ± 126 


340-363 


296 ± 22 


Xa(lP) 




568 


415 ± 170 


413-464 


386 ± 27 


h c (W) - 


>7] c (15)y 


909 


447 ± 272 




<600 


y{2S) 


Xco(lP)Y 


53.6 


21.4 ± 16.1 


26.0-40.3 


29.4 ± 1.3 


y(2S) -> 


Xcy{lP)Y 


45.2 


30.7 ± 13.6 


28.3-37.3 


28.0 ± 1.5 


y(2S) -> 


Xci{lP)Y 


31.6 


25.6 ± 9.5 


17.5-22.7 


26.5 ± 1.3 


r) c (2S) - 


>h c {\P)Y 


38.1 


31.0 ± 11.4 







Table 2: El decay rates for charmonium. Our pNRQCD results at LO, NLO (including error estimate) 
compared to a potential model calculation and the current PDG values [[24]]. 
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